Calculation of the Phase Field of a Vortex Pair 
on the Surface of a Multiplly Connected Substrate 
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The vortex pair phase field is calculated, through a solution of the Laplace equation, for a super- 
fluid film adsorbed on the surface of a model 3D connected porous medium. A number of different 
vortex-antivortex pair configurations are considered as functions of the porosity or an aspect ratio 
of the model porous structure. The calculated results give the energy of a pair for motion of the 
vortex cores around a pore or a solid rod as well as for the string configuration corresponding to a 
linear separation of the vortex pair through the porous medium. 
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The superfluid properties of thin 4 He films adsorbed 
on 3-D connected porous media have been studied for 
many years. These systems present a number of inter- 
esting physical phenomena: dilute bose gas behavior for 
low density 4 He adsorbed in Vycor glass [1]; 3-D crit- 
ical behavior observed for a wide range of systems [2]; 
and 2-D and 3-D crossover behavior as a function of sys- 
tem structure and pore size. Helium films adsorbed on 
large-scale flat substrates are known to exhibit the clas- 
sic 2-D Kosterlitz-Thouless phase transition [3-5] . In the 
case of 4 He films adsorbed on the 3-D connected sub- 
strate provided by a porous medium, it is clear that one 
is not dealing with ideal 2-D system; yet there remains 
the interesting question as to what extent vestigial as- 
pects of 2-D K-T behavior may still be present in the 
superfluid behavior of the system [6] . This question has 
attracted the attention of a number of theorists, the first 
being Minoguchi and Nagaoka (MN) [7] who argued for 
the existence of a 3-D superfluid phase transition for 3-D 
connected films and discussed vortex dynamics and the 
interrelation between 2-D and 3-D superfluid transitions. 
These authors pointed out that the energy of a vortex 
antivortex pair would have a string-like character with 
the energy growing linearly for large vortex separations. 
Subsequently, Machta and Guyer (MG) [8, 9] considered 
the 3-D porous media problem as well as the case of a 
film on a cylinder with infinite length. These authors, 
MN and MG, arrived at two main conclusions. First, 
that vortex pairs are confined by the string interaction 
and consequently a true vortex unbinding transition of 
the 2-D K-T type does not occur in a 3-D connected 
porous medium; and second, that the 3-D connectivity 
of the system enforces 3-D X-Y critical behavior at the 
superfluid transition. Gallet and Williams [10] were the 
first to emphasis that the string created by two vortices 
moving apart in the porous medium will not only create 
circulation around the solid surface over which the pair 
separates, but also around any loops through which a 
member of the pair moves. In all of the above papers the 
real complexity of the 3-D connected system is avoided 
through simplifying assumptions. 

In the present work we extend the work of the previ- 



FIG. 1: The porous substrate and the development of a rod. 
Every rod has the same dimensions. The lattice parameter 
and the width of the rod are defined in the figure, a = 2ro and 
a — nro/2 for comparisons with the cylinder model [9]. Then 
the porosity is 59%. The entire system has a translational 
symmetry in each three directions. 



ous authors with a more realistic calculation of the en- 
ergies for quantized vortex pairs taking into account the 
details of the 3-D porous structure. Following MN and 
MG, we employ the "Jungle Gym" model for the porous 
medium. The original model consisted of a cubic struc- 
ture composed of cylindrical rods with radius rg. In our 
case, we shall construct the lattice with rods of square 
cross section. 

Figure 1 shows a view of our lattice structure. An 
aspect ratio for this structure is given by, 7 = I /a where 
I is the lattice parameter (distance between nodes) and 
a is the width of the square rods. The ratio of the open 
pore volume to the total volume for a unit cell of the 
lattice defines the porosity, P = 1 — (37 — 2)/j 3 . The 
open volume of the pores forms a complementary cubic 
structure with the same lattice parameter and an aspect 
ratio, 7 = l/a = 7/(7 — 1), where a — I — a, is the 
length of a rod or the edge length of a pore. In making 
a comparison to the MG model, we shall initially choose 
a lattice parameter, I — (2 + 7r/2) ro and set a = nro/2. 
With this choice the length of the rods will be, a — 2ro 
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and the porosity has a value close to 0.59. 

In the present calculation the vortex pair phase field is 
computed on a 4x4x4 cubic lattice with 64 surface and 
interior lattice nodes and 192 rods. Periodic boundary 
conditions can be imposed that insure that the gradient 
in the phase normal to the surface of the cube is zero, 
i.e., the flow field induced by the vortex pair is restricted 
to the interior of the cube. A more ambitious project 
would be to perform these calculations on a larger lattice; 
for instance, one might take a 6x6x6 cubic structure 
with 216 nodes and 648 rods to be considered in the 
calculation. 

We shall assume that a uniform supcrfluid film with 
density per unit area, p, covers the surface of the lattice. 
In the present calculation, we will take into account the 
actual structure of the lattice in calculating the phase 
and superfluid flow fields for various configurations of a 
vortex-antivortex pair located on the surface of the lat- 
tice. For simplicity we start with the cores located close 
together on a single lattice rod. This rod will be cho- 
sen as close to the center of the lattice as possible given 
the discrete nature of the structure. The surface of the 
rod can be unfolded as illustrated in Figure 1. A coor- 
dinate r] indicates the distance between the vortex cores 
along the axis of the rod while £ gives the distance in the 
perpendicular direction (along the perimeter of the rod) . 
In the figure the vortex and antivortex cores are shown 
placed symmetrically about the center of the rod. The 
rod containing the vortex pair has a common boundary 
with four other adjacent rods and the phase field must 
be matched continuously across the common edges. 

At the start of the calculation, the initial position of 
the vortex pair is specified to be close together on an 
interior rod. The quantum phase, <f>, is calculated as a 
function of position on the surface of the lattice by a 
numerical solution of Laplace's equation on 64x64 grid 
on each rods while maintaining phase continuity across 
the matching edges of adjacent rods. Once the phase 
field is known, the superfluid velocity is obtained from 
the relation, v(r) = (K/mi)W4>(r) and the total energy 
from an integration of (l/2)pv 2 over the surface of the 
lattice. 

In Figure 2, three different configurations for a vor- 
tex pair are shown. In the first panel, 2a, the vortices 
are located on the rod labeled 1 and are separated by a 
distance, r\ = a/2. Phase values between and 2ir for 
this pair configuration are indicated by the color scheme. 
In panel 2b the color scale indicates the superfluid en- 
ergy distribution. When the vortex pair is restricted to 
a single rod there is very little flow on the adjacent rods. 
When the vortices are moved futher apart, until they are 
located at the corner intersection of the rods, the vortex 
flow field spills over onto the surfaces of the adjacent rods. 
In 2c the vortices are shown at the mid point of rods 2 
and 4. The vortices now face each other across the pore. 
As the vortices are moved apart, it is useful to keep track 
of the virtual core associated with the vortex-antivortex 
pair. The core will form a closed loop intersecting the 
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FIG. 2: A series of the calculation results. This solution is the 
example when vortex pair decouples along rj direction. The 
color corresponds to the phase from to 2n and it is periodical 
continuous in the panels a, c and d. (a) rj = 5/2. (b) The flow 
energy distribution of the configuration (a), (c) r\ = 22. The 
red broken line indicates the vortex core crossing the pore, 
(d) The symmetric flow pattern remaining after vortex pair 
annihilation. 



surface of the rods at the position of the vortices. The 
core then has two parts; one within the interior of the 
solid rods and another portion within the space of the 
open pores. When the separation is further increased to 
r] = 4a , the two elements of the pair will then overlap 
and the virtual core of the vortices will lie entirely within 
the rods. We are now left with a "core" vortex configura- 
tion, with a four-fold symmetric flow largely restricted to 
the circumference of the rods, 1 to 4. We shall call this 
flow configuration the rj - ring and describe the phase 
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FIG. 3: The vortex pair energy functions, 
(a) The energy function when a vortex pair expands along 
the 77 axis. The termination points are the annihilation points 
where a vortex pair annihilates, (b) The energy function when 
a vortex pair expands along the £ axis. The broken line is a 
normal K-T type logarithmic function, (c) The case when a 
vortex pair expand along the 77 direction through a pore [10]. 
The label A, B and C corresponds to the indicators in Fig. 5. 
The squares in the panels a and b are the 77 and £ ring energy 
respectively (see the main text). The broken lines in panels a 
and c are the analytic cylinder function [9] . 



configuration in Figure 2d. 

We now examine the change in the kinetic energy of 
the flow field as the pair separates along r\ direction. In 
Figure 3a we show a plot of the total energy as a function 
of the vortex-antivortex separation distance r\. In setting 
the energy scale, a supcrfluid mass per unit area is taken 
to be, p — 3.49 x 10~ 9 gm/cm 2 , so that a K-T transition 
causes at IK in the case of the infinitely spreaded 4 He 
film. The infinite cylinder calculation of MG[9] is shown 
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FIG. 4: The calculation results when the vortex pair expands 
in the £ direction, (a) £ = txtq. (b) The case of a pore vortex 
ring (see the main text. ). 



for comparison. It is interesting that these two calcula- 
tions agree rather well even as the vortex pair passes a 
corner in the lattice structure. 

We shall now consider pair separation in the £ direc- 
tion, perpendicular to 77, along the perimeter direction 
of the rod. The initial configuration is taken with the 
pair located at the center of the rod (77 = 0) separated 
by a short distance in the £ direction and the vortex pair 
moves on in the £ direction. When the separation is the 
perimeter distance, £ = 2nr , the pair will have encircled 
the rod and the cores will overlap and leave behind a 
continuous flow field along the axis of the rod. Figure 4 
shows the phase field on the surface of the lattice as the 
vortex pair expands in the £ direction. The development 
of the energy as a function of the separation parameter, 
£, is shown in panel b of Figure 3. The logarithmic de- 
pendence of the pair energy for a vortex-antivortex pair 
in 2-D is also shown as a function of pair separation. 
An interesting feature of the energy dependence is that 
a maximum occurs in the total energy shortly before the 
pair annihilates on the far side of the rod. After the an- 
nihilation, the system is left in a symmetric flow state 
similar to the flow pattern of a vortex ring in 3-D space; 
in the present case, however, the flow is restricted to the 
surface of the lattice. This remaining flow configuration 
had been called the "pore" vortex ring structure with 
the virtual vortex core threading the pores immediately 
surrounding the rod [8]. Here, we call this type of a pore 
vortex ring the £ - ring. 

A third configuration, shown in Figure 5, deserves con- 
sideration. The starting configuration is similar to the 
first case considered, with the vortices displaced in the 77 
direction. We shall take the initial position of the vor- 
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a nearly linear dependence of 7. 7=1 corresponds to the 
case that I = a; namely, there is no pore then no film let 
in the system. For comparison with this energy depen- 
dence it is worth considering a constant flow just on the 
surface of four adjacent rods forming a pore. We shall 
assume the magnitude of the velocity is fed by the phase 
difference of 2ir; i.e. v = h/Aam^. With this velocity, 
the total energy E=4x [{p/2)v 2 } x Aah = 12.6(7 - 1) 



(b) 

FIG. 5: A calculation result that vortex pair expands through 
pores. A vortex pair symmetrically expands from the middle 
point of an edge of the rod A. (a) The phase is described 
in color scheme. The indicators A, B and C correspond to 
the indicators in Figure 3. (b) The energy distribution of the 
same arrangement. At each corner where vortex have gone 
through, there are small energy concentration points. 

tices on the edge of the rod rather than in the middle 
of one face. Now, as 77 increases, the vortices can move 
in a straight line along the junction between rods and so 
on to the rods beyond. When a pair is separated with 
the elements moving past a junction in the lattice, from 
one rod to another, circulation is created not only around 
rods which have virtual core within but also about any 
loop in the lattice threaded by a single line of the core. 
Thus the energy of a separating vortex pair does not only 
grow linearly with separation in the string fashion but 
increases in a step-like fashion whenever a vortex moves 
through a pore. The result is an even stronger "confine- 
ment" [7] than results from the simple string interaction. 

We have also calculated the ring energies for different 
values of 7 (or porosoity) . The results are shown in Fig- 
ure 6. ^ ring is the £ - ring energy while E^ ing the 77 - ring 
energy. From symmetrical point of view, E^ infl for 30% 
porosity is E^ ing for 70% porosity and vice versa; in other 

words E^ ing for a certain 7 is E^ ing for the same 7 value 
as 7. For small values of 7, the r\ - ring energy is nearly a 
logarithmic function of 7, while for large 7, it approaches 



FIG. 6: The rj and £ ring energies are shown as open and 
closed circles, respectively, as functions of the aspect ratio 7. 



(K), where we have again set p = 3.49 x 10 -9 gm/cm 2 
in setting the energy scale. This energy function is a 
little smaller than E^ ing but parallel for large values of 
7. This difference is caused by the fact that the actual 
velocity distribution is not constant because of the 3D 
connectivity. 

The energy dependence tells us whether a pore vortex 
ring could be excited or not. Around 7 = 2, both E^ irag 

and E^ ing are too high to be excited and a pore vor- 
tex ring would be rarely excited. Because we now know 
the energy behaves like a cylinder case in both £ and 
■q direction until it approaches the maximum, we could 
constract the renomalization calculation by an imitation 
energy function, which expresses a vortex string energy 
for both axes. The study would be closely connected to 
the experimental study on the Vycor glass substrate [2] 
as well as on the other Porous glass substrates [6] . As 7 
increases, the porosity becomes higher and a pore vortex 
ring comes to be excited more easily. Then some sim- 
ple mean field theory would better describe the physical 
property of the system like 4 He film on aero gel [11]. 

The authors thank many people. We thank Prof. J.D. 
Reppy for fruitful discussions and suggestions. The col- 
ored graph was described using the basic graphics class 
library for C by Mr. K. Yamada[12]. 
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